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Live-fire  military  training  can  deposit  millimeter-sized  particles  of  high  explosives  (HE)  on  surface  soils 
when  rounds  do  not  explode  as  intended.  Rainfall-driven  dissolution  of  the  particles  then  begins  a  pro¬ 
cess  whereby  aqueous  HE  solutions  can  enter  the  soil  and  groundwater  as  contaminants.  We  dripped 
water  onto  individual  particles  of  TNT,  Tritonal,  Comp  B  and  Octol  to  simulate  how  surface-deposited 
HE  particles  might  dissolve  under  the  action  of  rainfall  and  to  use  the  data  to  verify  a  model  that  predicts 
HE  dissolution  as  a  function  of  particle  size,  particle  composition  and  rainfall  rate.  Particle  masses  ranged 
from  1.1  to  17  mg  and  drip  rates  corresponded  to  nominal  rainfall  rates  of  6  and  12  mm  h-1.  For  the  TNT 
and  Tritonal  particles,  TNT  solubility  governed  dissolution  time  scales,  whereas  the  lower-solubility  of 
RDX  controlled  the  dissolution  time  of  both  RDX  and  TNT  in  Comp  B.  The  large,  low-solubility  crystals 
of  HMX  slowed  but  did  not  control  the  dissolution  of  TNT  in  Octol.  Predictions  from  a  drop-impingement 
dissolution  model  agree  well  with  dissolved-mass  timeseries  for  TNT,  Tritonal  and  Comp  B,  providing 
some  confidence  that  the  model  will  also  work  well  when  applied  to  the  rainfall-driven,  outdoor  disso¬ 
lution  of  these  HE  particles. 

Published  by  Elsevier  Ltd. 


1.  Introduction 

Military  ranges  provide  soldiers  the  opportunity  to  train  using  a 
variety  of  munitions.  However,  live-fire  training  can  deposit  high 
explosives  (HE)  such  as  2,4,6-trinitrotoluene  (TNT),  1,3,5-hexahy- 
dro-l,3,5-trinitrotriazine  (RDX)  and  octahydro-l,3,5,7-tetranitro- 
1,3,5,7-tetrazocine  (HMX)  on  the  range.  Detonations  scatter  HE 
particles  broadly  over  surface  soils.  High-order  detonations  scatter 
pm-size  HE  particles  and  low-order  (LO)  detonations  scatter  mm- 
to-cm-sized  HE  particles  (Taylor  et  al.,  2006).  Given  the  low  drink¬ 
ing  water  screening  levels  for  TNT  (2.2  pgL-1)  and  RDX 
(0.6  pgL-1)  (  EPA,  2008)  gram-to-kilogram  quantities  of  these 
explosives  in  aqueous  solution  may  contaminate  an  aquifer,  as  oc¬ 
curred  at  the  Massachusetts  Military  Reservation  (Clausen  et  al., 
2006).  The  low  drinking  water  levels  reflect  their  toxicity  to  human 
health  (ATSDR,  1995). 

Range  characterization  studies  summarized  by  Pennington  et  al. 
(2006)  found  HE  in  surface  soils  at  all  27  military  installations  sam¬ 
pled  and  HE  was  found  in  groundwater  at  Fort  Lewis  WA  (Jenkins 
et  al.,  2001),  and  at  three  Canadian  installations  (Martel  et  al., 
1999;  Pennington  et  al.,  2006).  Dissolution  of  HE  particles  by  pre¬ 
cipitation  has  to  occur  before  aqueous  transport  or  bio-degradation 
of  the  HE.  Therefore  dissolution  rates  of  HE  particles  as  a  function 
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of  particle  size,  composition  and  rainfall  rate  are  needed  to  predict 
HE  aqueous  influx  to  surface  soils  at  training  ranges. 

Lynch  et  al.  (2002a, b)  and  Phelan  et  al.  (2003)  conducted, 
respectively,  stirred-bath  and  glass-bead-column  experiments  to 
measure  the  aqueous  dissolution  rates  of  collections  of  HE  parti¬ 
cles.  We  followed  the  method  described  in  Lever  et  al.  (2005) 
and  dripped  water  on  individual,  field-collected,  mm-sized  parti¬ 
cles  of  four  common  high-explosive  formulations:  TNT,  Tritonal 
(80:20  mix  of  TNT  and  aluminum  flakes),  Composition  B  (59:39:1 
mix  of  RDX:TNT:wax)  and  Octol  (70:30  mix  of  HMX:TNT).  These 
are  all  melt  cast  explosives  where  aluminum,  RDX  and  HMX  crys¬ 
tals  are  added  to  molten  TNT.  The  experiments  simulated  the  rain- 
fall-driven  dissolution  of  HE  particles  scattered  on  surface  soils.  We 
then  used  the  resulting  dissolved-mass  timeseries  to  validate  a 
drop-impingement  dissolution  model  developed  for  Comp  B  (Lever 
et  al.,  2005).  By  expanding  the  types  of  HE  formulations  tested  and 
using  two  simulated  rainfall  rates,  our  intention  was  to  improve 
confidence  in  the  model’s  predictions  of  rainfall-driven,  dissolution 
of  HE  particles  on  training  ranges. 


2.  Experimental  methods 

2 A.  HE  Particle  samples 

TNT  is  a  single-component  explosive,  as  is  Tritonal  albeit  with 
20%  aluminum  flakes,  while  Comp  B  and  Octol  are  two-component 
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explosives.  We  tested  two  particles  for  each  HE  formulation,  all 
collected  from  training  ranges.  The  TNT  samples  were  collected 
after  a  LO  detonation  of  a  155-mm  round,  the  Comp  B  samples 
were  from  the  LO  detonation  of  an  81 -mm  round,  and  the  Tritonal 
samples  were  from  a  LO  detonation  of  a  230-kg  (500-lb)  bomb.  We 
obtained  a  single  Octol  particle  from  soils  sampled  on  an  anti-tank 
firing  range  and  split  it  to  form  the  two  samples  tested  here. 

Before  starting  the  drip  experiments,  we  individually  weighed 
the  dry  HE  particles  on  a  Mettler  Toledo  MX5  microbalance  (uncer¬ 
tainty  ±  0.01  mg)  and  photographed  them  using  a  camera  attached 
to  a  binocular  microscope.  Part  way  through  the  tests,  we  re-pho¬ 
tographed  the  particles  to  document  changes  to  their  appearance. 

Lever  et  al.  (2005)  showed  that  mm-sized  Comp  B  detonation 
residues  varied  substantially  in  their  RDX/TNT  ratios,  averaging 
1.74  ±  0.28  (±1  standard  deviation)  compared  with  a  nominal  com¬ 
position  of  60/39  =  1.54.  Significant  deviations  in  composition  are 
likely  for  field-recovered  HE  particles  with  component  crystals  that 
are  large  compared  with  overall  particle  size.  To  measure  this  var¬ 
iability  for  Tritonal,  we  weighed  21  mm-sized  particles  on  a  Sarto- 
rius  electronic  balance  (uncertainty  ±  0.1  mg),  dissolved  them  in 
acetonitrile  and  analyzed  their  TNT  compositions.  We  did  not  have 
sufficient  Octol  to  investigate  variability  in  Octol  particle 
compositions. 

2.2.  Laboratory  setup 

Two  Cole  Palmer  multi-syringe  pumps  allowed  us  to  drip  de¬ 
ionized  water  simultaneously  onto  eight  separate  particles,  two 
particles  for  each  of  the  four  HE  formulations  at  two  different  drip 
rates.  A  needle  on  each  syringe  dripped  water  onto  an  HE  particle 
resting  on  a  10-mm-diameter  glass  frit  at  the  base  of  a  Buchner 
funnel  that  drained  into  a  20-  or  40-mL  scintillation  vial.  The  vials 
were  replaced  daily  and  the  volume  and  HE  concentration  of  the 
effluent  analyzed.  These  data  provided  dissolved  HE  mass  as  a 
function  of  time  for  each  of  the  eight  particles.  At  the  end  of  the 
tests,  we  extracted  each  frit  and  funnel  with  acetonitrile  and  ana¬ 
lyzed  the  solutions  for  HE  mass. 

The  programmed  flow  rates  for  each  pump  were  0.5  and 
1.0  mLh-1.  We  measured  the  actual  flow  rates  and  number  of 
drops  per  hour  for  each  syringe.  These  averaged  0.47  ± 
0.06  mL  h-1  with  26  drops  h-1  and  0.95  ±  0.09  mL  h_1  with  56 
drops  h-1  at  the  two,  programmed  rates.  For  reference,  these  flow 
rates  correspond  to  rainfall  rates  of  5.9  and  12.0  mm  h_1  based  on 
the  cross-sectional  area  of  each  glass  frit.  The  average  volume  of 
the  arriving  drops  was  larger  than  that  expected  for  the  natural 
rainfall  at  the  same  rates  (Pruppacher  and  Klett,  1997)  and  each 
arriving  drop  wetted  the  particle  on  the  frit.  The  de-ionized  water 
used  for  the  tests  was  at  room  temperature  and  the  latter  was 
fairly  constant  throughout  the  tests  at  22  ±  1  °C.  The  tests  were 
normally  halted  over  the  weekends  and  the  particles  allowed  to 
dry. 


2.3.  Analytical  methods 

Explosive  concentrations  were  determined  following  SW-846 
Method  8330B  (EPA,  2006).  One  mL  of  the  effluent  was  added  to 
2  mL  of  de-ionized  water  and  to  1  mL  of  acetonitrile  and  then  fil¬ 
tered  through  a  0.45  pm  Millipore  syringe  filter.  High  Performance 
Liquid  Chromatography  separated  TNT,  RDX,  HMX,  and  their  co¬ 
contaminants  and  breakdown  products.  We  used  a  Waters  NovaP- 
ak  C8  column  (3.9  mm  by  150  mm)  eluted  at  1.4  mL  min-1  (28  °C) 
with  85:15  water: isopropanol  mix  and  detected  by  UV  at  254  nm. 
A  commercially  available  standard  (Restek)  was  used  for  calibra¬ 
tion.  We  prepared  1-ppm  and  10-ppm  dilutions  of  these  standards. 
The  1  ppm  standards  were  run  every  ten  samples  and  blanks  were 
run  before  and  after  each  standard  run.  The  10  ppm  standards 
were  interspersed  with  the  samples  as  internal  checks,  with  blanks 
after  each  to  prevent  carryover.  Based  on  the  concentrations  of  the 
standards  and  the  precision  of  the  effluent  volumes,  we  estimate 
that  the  cumulative  dissolved  HE  masses  have  uncertainties  of 
about  ±1%. 


3.  Results 

3.1.  He  mass  recovery 

We  analyzed  21  Tritonal  particles  for  compositional  variability. 
The  percentage  of  TNT  by  weight  ranged  from  52%  to  84%  of  parti¬ 
cle  mass  and  averaged  77.5%  ±  6.6%.  Thus,  the  average  of  the  21 
particles  fell  well  within  1  standard  deviation  of  the  nominal  com¬ 
position  of  Tritonal  (80%  TNT). 

Table  1  lists  the  dissolution  test  parameters  and  the  recovered 
mass  for  the  eight  drip  test  particles.  The  nominal  drip  rates  were 
0.5  mL  h-1  and  1 .0  mL  h_1  for  the  “1”  and  “2”  particle  designations, 
respectively.  We  recovered  essentially  100%  of  the  mass  for  the 
two  TNT  particles  and  for  Comp  B  2.  The  recovered  HE  mass  for 
Comp  B  1  (96.5%)  was  close  to  the  average  (96%  ±1%)  found  by  Le¬ 
ver  et  al.  (2005)  for  30  individual  Comp  B  particles  obtained  from 
the  same  low-order  detonation.  The  two  particles  dissolved  here 
had  RDX/TNT  ratios  of  1.69  and  1.87,  similar  to  the  average  ratio 
obtained  by  Lever  et  al.  (2005)  of  1.74  ±0.28.  We  measured  less 
than  1%  HMX,  an  impurity  in  RDX,  in  these  Comp  B  particles. 

The  recovered  TNT  mass  was  67.8%  and  78.5%  of  the  initial  mass 
for  Tritonal  1  and  Tritonal  2,  respectively.  This  included  TNT  ex¬ 
tracted  from  the  frits.  Both  particles  fell  within  2a  of  the  nominal 
Tritonal  composition.  Aluminum  grains  constituted  the  remaining 
mass  for  the  Tritonal  particles. 

For  Octol  1  and  Octol  2,  we  recovered  1.86  mg  and  3.61  mg  of 
TNT,  representing  28.4%  and  20.8%  of  the  initial  particle  masses, 
respectively.  Octol  nominally  contains  30%  TNT,  but  we  would  ex¬ 
pect  millimeter-sized  particles  to  vary  significantly  from  this  va¬ 
lue  owing  to  the  large  (~1  mm)  crystals  of  HMX  present  in  the 
TNT  matrix.  Both  dissolution  tests  continued  for  many  days  with 


Table  1 

Dissolution  test  parameters  and  the  recovered  mass  for  HE  particles. 


All  the  TNT  mass. 


1076 


S.  Taylor  et  al.  f  Chemosphere  75  (2009)  1074-1081 


TNT  1 :  5.34  mg 


TNT  2:  9.59  mg 


Trit  1-  1.89  mg 


Trit  2-  6.40  mg 


Comp  B  1-  2.31  mg 


Comp  B  2-  9.09  mg 


Octol  I  -  6.53  mg 


Octol  2-  17.33  mg 


TNT  1  after  30  days 


TNT  2  after  21  days 


1  mm 

Trit  1  after  30  days 


Comp  B  I  after  30  days 


□1 


1  mm 


Comp  B  2  after  2 1  days 


I  mm 

Octol  I  after  30  days 


1,  mm 

Octol  2  after  2 1  days 


1  mm 


Fig.  1.  Photos  and  initial  masses  of  mm-sized  particles  used  in  the  laboratory  tests 
paired  with  photos  of  the  same  particle  taken  partway  through  the  tests. 


no  measurable  TNT  in  the  water  samples,  and  we  found  no  addi¬ 
tional  TNT  in  a  series  of  acetonitrile  extractions  of  the  funnels  and 
frits  at  the  end  of  the  tests.  Thus,  we  think  that  the  dissolution 


tests  recovered  all  the  TNT  initially  present  in  the  two  Octol 
particles. 


3.2.  Appearance  of  test  particles 

Fig.  1  shows  photographs  of  the  8  particles  taken  at  the  begin¬ 
ning  and  partway  through  the  dissolution  tests,  at  30  days  for  the 
0.5  mLh-1  tests  and  at  21  days  for  the  1.0  mLh-1  tests.  The  TNT 
particles  became  smoother  and  smaller  but  retained  their  original 
shapes.  The  Tritonal  particles  became  smaller  and  slightly  bumpier 
as  TNT  dissolved  exposing  the  aluminum  grains.  The  Comp  B  par¬ 
ticles,  on  the  hand,  became  noticeably  bumpier  and  ‘sugary’  look¬ 
ing  as  dissolution  of  the  surface  TNT  revealed  the  larger 
(~  0.1  mm),  slower-dissolving  RDX  crystals.  The  Octol  particles 
disaggregated  completely  when  their  TNT  matrix  dissolved  leaving 
large  (~1  mm)  HMX  crystals  (see  Fig.  1,  Octol  2  at  21  d). 


3.3.  Dissolution  modeling  and  comparison  with  dissolved-mass 
timeseries 

Lever  et  al.  (2005)  developed  and  validated  a  drop-impinge¬ 
ment  model  for  the  rainfall-driven  dissolution  of  mm-sized  Comp 
B  particles.  It  pertains  to  the  practical  case  where  spatially  isolated 
HE  particles  reside  on  well-draining  surface  soils  and  thus  are  ex¬ 
posed  to  direct  impingement  by  raindrops.  The  present  experi¬ 
ments  simulated  these  circumstances.  With  rare  exceptions, 
water  quickly  drained  through  the  glass  frits  so  that  impinging 
drops  interacted  directly  with  the  HE  particles.  The  present  exper¬ 
iments  offer  additional  data  to  validate  this  model  for  other  HE  for¬ 
mulations  and  for  two  simulated  rainfall  rates. 

The  drop-impingement  dissolution  model  assumes  that  rain¬ 
drops  hit  and  wet  an  HE  particle.  Between  raindrops,  the  particle 
holds  a  stagnant  water-layer  against  its  surface,  which  saturates, 
via  diffusion,  with  HE.  Arrival  of  the  next  raindrop  washes  away 
the  dissolved  HE  and  refreshes  the  stagnant  layer.  The  dissolution 
rate,  m,  (g  s_1)  of  HE  species  j,  averaged  over  a  drop-arrival  interval, 
td  (s),  is  thus 


where  S,  is  the  solubility  of  species  j  in  water  (g  cm-3)  and  V\  is  the 
water-layer  volume  (cm3).  For  a  spherical  particle,  a  water-layer  of 
thickness  h  (cm)  will  effectively  saturate  (average  concentra¬ 
tion  >  0.9  Sj)  provided 

td>^  =  ts  (2) 

where  Dj  is  the  diffusion  coefficient  of  species  j  into  water  (cm2  s-1) 
and  ts  defines  the  layer  saturation  time.  Additionally,  for  Eq.  (1)  to 
apply,  the  volume  of  an  arriving  water  drop,  Vd,  must  be  larger  than 
the  volume  of  the  stagnant  water-layer,  which  generally  limits  its 
applicability  to  mm-sized  and  smaller  particles. 

Each  HE  species  in  a  particle  could  dissolve  independently,  fol¬ 
lowing  Eq.  (1 ),  if  it  is  in  contact  with  the  water-layer  and  has  a  suf¬ 
ficiently  high  diffusion  coefficient  to  satisfy  Eq.  (2).  However,  Comp 
B  and  Octol  particles  consist  of  low-solubility  crystals  (RDX  and 
HMX,  respectively)  embedded  in  higher-solubility  TNT  matrixes. 
Initially,  TNT  on  the  surface  of  these  particles  can  dissolve  indepen¬ 
dently  into  the  water-layer,  leading  to  high  initial  TNT  dissolution 
rates.  However,  the  lower-solubility  crystals  eventually  restrict  the 
pathways  for  TNT  molecules  to  diffuse  into  the  water-layer.  TNT 
diffusion  into  the  layer  then  becomes  analogous  to  molecular  dif¬ 
fusion  through  a  porous  medium,  where  porosity  and  tortuosity 
combine  to  reduce  significantly  the  effective  diffusion  coefficient 
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Table  2 

Model  parameters.  Particle  densities  are  from  Army  Materiel  Command  (1971 ).  Aqueous  solubilities  are  from  Lynch  et  al.  (2001 )  at  22  °C,  and  diffusion  coefficients  are  from  Lynch 
et  al.  (2002b)  at  25  °C;  listed  are  TNT  values  for  TNT  and  Tritonal,  RDX  values  for  Comp  B,  and  HMX  values  for  Octol.  The  two  densities  for  Tritonal  reflect  the  different  TNT / 
aluminum  ratios  of  the  two  particles. 


0  20  40  60  80  100  120 


Time  (d) 

Fig.  2.  Dissolved-mass  versus  time  for  the  TNT  test  particles  (symbols  are  measured  values,  smooth  curves  are  modeled  values). 


(Bear,  1972;  Chambre  and  Pigford,  1984).  Consequently,  TNT  will 
not  saturate  the  water-layer  between  drops  and  its  effective  disso¬ 
lution  rate  will  decrease. 

Lever  et  al.  (2005)  observed  this  effect  during  laboratory  drip 
tests  of  Comp  B  particles.  They  implemented  the  limiting  case 
where  RDX  dissolution  controls  the  dissolution  of  TNT  by  using 
Eq.  (1)  for  RDX  and  assuming  that  TNT  dissolves  at  a  rate  that 
maintains  the  bulk  RDX/TNT  mass  ratio: 


miNT  =  fllRDX  TNT  (3) 

Prdx 

where  pTnt  and  pRDX  are  the  mass  densities  of  each  species  in  the 
original  Comp  B  particle.  This  approach  ensured  that  the  time  scale 
for  RDX  dissolution  controlled  the  predicted  TNT  dissolution  rate. 
We  applied  this  approach  for  both  Comp  B  and  Octol  particles, 
substituting  pHMx  for  pRDX  in  Eq.  (3)  for  Octol. 
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We  implemented  the  drop-impingement  model  for  the  eight 
test  particles  by  treating  each  particle  as  a  sphere  of  equivalent 
mass,  so  that  the  water-layer  volume  is  simply 

Vl  =  *n[(a  +  h)3  -a3}  (4) 

where  a  is  sphere  radius  (cm).  Because  we  could  not  measure  the 
water-layer  thicknesses,  we  determined  h  by  fitting  the  predictions 
to  the  measured  dissolved-mass  timeseries  for  the  slower-dissolv¬ 
ing  species  in  each  particle.  Effects  of  non-spherical  particle  shape 
are  thereby  folded  into  h.  Table  2  summarizes  the  parameters 
needed  to  run  the  drop-impingement  model  and  to  check  its  valid¬ 
ity  (td  >  ts.  Vd  >  Vi). 

Figs.  2  and  3  show  the  measured  and  predicted  dissolved- 
masses  for  the  TNT  and  Tritonal  particles,  respectively.  These  par¬ 
ticles  contain  only  TNT,  so  Eq.  (1)  applies  directly.  The  model  pre¬ 
dicts  the  dissolution  of  all  four  particles  quite  well.  The  average 
layer  thicknesses  for  the  TNT  and  Tritonal  particles  were  essen¬ 
tially  the  same  at  0.085  ±  0.014  mm  and  0.086  ±  0.006  mm,  respec¬ 
tively.  Note  that  for  the  particles,  the  drop-arrival  times  were 
longer  than  the  water-layer  saturation  times,  and  the  drop  vol¬ 
umes  were  larger  than  the  water-layer  volumes  (Table  2)  as  re¬ 
quired  by  Eq.  (1). 


Fig.  4  shows  the  measured  and  predicted  dissolved-masses  for 
the  two  Comp  B  particles.  The  predicted  RDX  dissolution  agrees 
quite  well  with  the  measured  results  for  both  particles.  In  both 
tests,  TNT  initially  dissolved  quickly  but  eventually  was  controlled 
by  the  slower  dissolution  of  RDX.  The  model  predicts  this  behavior 
reasonably  well  using  the  limiting  case  of  RDX  control  of  TNT  dis¬ 
solution  (Eq.  (3)).  By  comparison,  it  predicts  much  faster  dissolu¬ 
tion  under  the  assumption  that  TNT  dissolves  independently  of 
RDX.  Interestingly,  the  fastest  dissolution  rates  observed  near  the 
beginning  of  the  tests  do  approximate  those  for  independent  TNT 
dissolution.  The  transition  from  independent  TNT  dissolution  to 
RDX-controlled  dissolution  probably  reflects  restriction  of  TNT  dif¬ 
fusion  into  the  water-layer  by  the  porous  medium  of  RDX  crystals. 

The  fitted  water-layer  thicknesses  average  0.101  ±  0.013  mm 
for  the  two  Comp  B  particles  (Table  2).  This  compares  well  with 
the  average  value  (0.13  ±  0.05  mm)  obtained  for  the  four  Comp  B 
particles  tested  by  Lever  et  al.  (2005).  For  both  Comp  B  particles, 
the  drop  volumes  were  larger  than  the  water-layer  volumes,  and 
the  drop-arrival  times  were  longer  than  the  RDX  saturation  times 
(Table  2).  TNT  saturation  times  must  increase  to  account  for  even¬ 
tual  RDX  control  of  TNT  dissolution. 

Fig.  5  shows  the  measured  and  predicted  dissolved-masses  for 
the  two  Octol  particles.  As  with  Comp  B,  measured  TNT  dissolution 


Fig.  3.  Dissolved-mass  versus  time  for  the  Tritonal  test  particles  (symbols  are  measured  values,  smooth  curves  are  modeled  values). 
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Fig.  4.  Dissolved-mass  versus  time  for  the  Comp  B  test  particles  (symbols  are  measured  values,  smooth  curves  are  modeled  values).  Dissolution  of  RDX  eventually  controls 
dissolution  of  TNT  (modeled  via  Eq.  (3)).  Independent  TNT  dissolution  would  occur  much  faster. 


initially  proceeded  quickly,  with  maximum  dissolution  rates 
approaching  the  predicted  rates  for  independent  TNT  dissolution. 
Unlike  Comp  B,  however,  the  larger  HMX  crystals  did  not  eventu¬ 
ally  control  TNT  dissolution,  and  all  the  TNT  dissolved  from  both 
particles,  leaving  HMX  crystals.  We  terminated  both  tests  after 
many  days  with  no  TNT  in  the  water  samples.  The  fitted  water- 
layer  thicknesses  averaged  0.155  ±  0.049  mm,  much  larger  than 
those  for  the  other  particles  (Table  2).  Nevertheless,  the  calculated 
saturation  times  (Eq.  (2))  based  on  either  TNT  or  HMX  diffusion 
were  less  than  the  drop-arrival  times,  and  the  water-film  volumes 
were  less  than  the  drop  volumes,  as  required  by  the  model. 

4.  Discussion 

The  laboratory  tests  dripped  water  on  individual  HE  particles 
collected  from  field  detonations.  They  more  closely  simulate  the 
physical  circumstances  of  HE  particles  exposed  to  rainfall  on  sur¬ 
face  soils  than  do  stirred-bath  (Lynch  et  al.,  2002a, b)  or  glass-bead 
column  experiments  (Phelan  et  al.,  2003).  They  yield  estimates  of 
aqueous  HE  influx  to  surface  soils  as  functions  of  particle  size,  par¬ 
ticle  composition  and  rainfall  rate  while  avoiding  the  complexity  of 
aqueous-phase  HE-soil  interactions  that  occur  during  soil-column 
experiments  (Pennington  et  al.,  2006;  Morley  et  al.,  2006).  TNT 


and  RDX  mass  recoveries  from  all  particles  were  very  good;  HMX 
mass  recovery  was  complicated  by  its  low-solubility  and  disaggre¬ 
gation  of  the  Octol  particles  following  complete  TNT  dissolution. 

Changes  in  the  appearance  of  the  particles  provided  insight  into 
their  dissolution  mechanics.  The  single  HE-species  particles,  TNT 
and  Tritonal,  essentially  retained  their  original  morphologies  as 
they  dissolved,  and  their  dissolution  timeseries  are  fairly  simple 
curves.  The  concentration  of  aluminum  grains  in  Tritonal  is  appar¬ 
ently  too  low  to  disrupt  dissolution  of  the  TNT.  By  comparison,  the 
Comp  B  and  Octol  particles  quickly  became  lumpy  as  surface  TNT 
dissolved  to  reveal  underlying  RDX  and  HMX  crystals.  Subse¬ 
quently,  the  Comp  B  particles  became  smaller  conglomerations  of 
RDX  crystals  shielding  internal  TNT,  while  the  Octol  particles  lost 
all  their  TNT  and  disaggregated  into  collections  of  HMX  crystals. 

We  do  not  understand  the  reasons  for  discontinuities  in  the 
measured  dissolution  rates  for  some  particles  (e.g.,  TNT  1,  Tritonal 
2  and  Comp  B  1 ).  The  discontinuities  do  not  coincide  with  the  dates 
of  particle  removal  for  photographing  or  observed  splitting  of  the 
particles,  the  good  mass  balances  confirm  that  the  particles  did 
not  lose  pieces  during  handling,  and  several  particles  did  not  dis¬ 
play  discontinuities  beyond  random  fluctuations. 

The  drop-impingement  dissolution  model  predicts  very  well  the 
dissolved-mass  timeseries  of  the  TNT  and  Tritonal  particles.  The 
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Fig.  5.  Dissolved-mass  versus  time  for  the  Octol  test  particles  (symbols  are  measured  values,  smooth  curves  are  modeled  values).  TNT  dissolution  rate  apparently  was  neither 
controlled  by,  nor  independent  of,  HMX  dissolution. 


single  tuning  parameter,  water-layer  thickness,  h ,  was  reasonably 
consistent  among  the  four  particles,  averaging  0.086  ±  0.009  mm. 
The  model  also  predicts  well  the  overall  dissolution  behavior  of 
the  Comp  B  particles  and  provides  insight  into  the  governing  pro¬ 
cesses.  Initially  TNT  dissolves  independently  of  RDX,  but  increasing 
restriction  of  TNT  diffusion  by  the  RDX  crystals  eventually  causes 
TNT  dissolution  to  be  controlled  by  RDX  dissolution.  With  addi¬ 
tional  work,  we  hope  to  quantify  this  effect  by  modeling  TNT  diffu¬ 
sion  within  a  changing  porous  RDX  medium. 

The  model  does  not  predict  well  the  TNT  dissolved-mass 
timeseries  for  the  Octol  particles.  The  large,  low-solubility  HMX 
crystals  appeared  to  impede  but  not  control  TNT  dissolution. 
That  is,  neither  limiting  model  case  (independent  or  HMX-con- 
trolled  dissolution)  fits  the  observed  TNT  data.  The  reduced 
influence  of  HMX  on  TNT  dissolution,  compared  with  the  con¬ 
trolling  role  of  RDX  in  Comp  B,  is  consistent  with  less  constraint 
of  TNT  diffusion  into  the  surrounding  water-layer.  The  higher  fit¬ 
ted  values  of  layer  thickness  (0.155  ±  0.049  mm)  could  reflect  the 
creation  of  surface  area  as  the  Octol  particles  disaggregated  fol¬ 
lowing  TNT  loss.  These  values  give  reasonable  agreement  to  the 
limited  HMX  dissolution  data,  but  more  work  is  needed  to  quan¬ 
tify  the  impeding  effect  of  HMX  crystals  on  TNT  dissolution  in 
Octol. 


The  drop-impingement  model  accounts  for  differences  in  rain¬ 
fall  rates  via  differences  in  average  drop-arrival  time,  td,  in  Eq. 
(1).  The  generally  good  agreement  between  the  observed  and  pre¬ 
dicted  dissolution  data  at  two  simulated  rainfall  rates  suggests  that 
td  does  indeed  scale  dissolution  rate. 

We  have  initiated  a  series  of  outdoor  dissolution  tests  where 
cm-sized  chunks  of  HE  are  exposed  to  natural  weather  conditions. 
These  tests  should  closely  simulate  the  dissolution  of  isolated  HE 
pieces  on  range  soils  and  provide  a  more  realistic  validation  of 
the  drop-impingement  model.  We  are  developing  a  “large-particle” 
version  of  the  model  to  circumvent  the  restriction  that  raindrop 
volume  be  greater  than  water-film  volume.  Nevertheless,  the  mod¬ 
el  predicts  quite  well  the  dissolution  of  millimeter-size  TNT,  Tri- 
tonal  and  Comp  B  particles.  This  provides  confidence  it  will  also 
work  well  when  applied  to  the  rainfall-driven,  outdoor  dissolution 
of  these  HE  materials. 
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